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Abstract: The present study revisited the controversial taxonomic status of Petaurista yunanensis, P. philippensis, P. 
hainana, and P. petaurista by using a considerably extended set of morphometrical characters (26 cranial variables from 
60 adult specimen skulls). The results revealed no sexual dimorphism in any of the four species but confirmed significant 
craniometric differences among the four species in both the principal components analysis (PCA) and discriminant 
function analysis (DFA), with the greatest distinction observed between P. petaurista and other Petaurista species. Both 
univariate and multivariate analysis indicated that the morphological differences between P. yunanensis and P. 
philippensis were less than that between P. philippensis and P. hainana. The morphometric results were concordant in 
geographic patterns with mtDNA data from previous studies and indicated that P. petaurista, P. hainana, P. philippensis, 
and P. yunanensis could be recognized as valid species. 
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Giant flying squirrels (Petaurista) occur in lowlands Because various species and subspecies with significant 
and mountains up to 4 000 m in East and Southeast Asia geographical variations are included within this genus, 
from Pakistan and Kashmir. Currently, eight to ten forms the taxonomy and the intra- and inter-specific 
are commonly accepted as valid Petaurista species, with phylogenetic relationships remain unclear and 
each divided into various forms or subspecies (Corbet & inconclusive. Several taxonomic studies and more than 
Hill, 1992; Wang, 2003; Thorington & Hoffmann, 2005). 18 Petaurista forms, subspecies, or species have been 
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described on the basis of dental and cranial 
characteristics and external structures (Allen, 1940; 
Corbet & Hill, 1992; Ellerman, 1940; Hoffmann et al, 
1993; Wang, 2003; Zhang et al, 1997). 

Petaurista are widely distributed in China and more 
than ten distinct species are recognized (Corbet & Hill, 
1992; Wang, 2003; Zhang et al, 1997). However, several 
of these species are referenced with very few specimens 
or based solely on skins with no corresponding skulls 
(Allen, 1940; Ellerman, 1940), and some are actually the 
synonyms or subspecies of either the P. petaurista 
complex or the P. philippensis complex due to intraspecific 
geographic variations across their distributions in Asia. 

Corbet & Hill (1992) treated P. albiventer in 
Pakistan and southwest China as the synonym of P. 
petaurista and recognized P. philippensis as a distinct 
species consisting of many forms formerly assigned to P. 
petaurista, including forms distributed in Taiwan (P. 
grandis), southwest Yunnan (P. yunanensis), and Hainan 
(P. hainana). After comparing the pelage and cranial 
characteristics of P. petaurista and P. hainana, Huang et 
al (1995) considered P. hainana to be a valid species, but 
Wang (2003) treated P. hainana as a subspecies of P. 
yunanensis. Thorington & Hoffmann (2005) treated all 
Petaurista forms as eight valid species instead of nine as 
suggested by Corbet & Hill (1992), but they accepted the 
specific validity of P. philippensis and the subspecies 
status of P. yunanensis and P. hainana. Patterns of 
genetic variations observed in the complex of P. 
philippensis based on cytochrome b genes indicated that 
P. hainana, P. albiventer, and P. yunanensis could be 
distinct species (Yu et al, 2006). Some forms included in 
P. philippensis warranted separate specific rank based on 
molecular data (Oshida et al, 2000a, b; Yu et al, 2006), 
but without further evidence from morphometric data, 
much remains to be done to ascertain conclusively these 
specific conclusions. 

Most recent phylogenetic studies have focused on 
molecular data analysis, but tracing changes in 
morphological characters is also an important way to 
evaluate the distribution of the characters on which those 
taxonomic units are based. Morphometric data are 
important to understand biological phenomena and have 
dental, and body 
measurements of many mammals (Mufioz-Mufioz & 
Perpinan, 2010; Slabova & Frynta, 2007; Zelditch et al, 
2004). Quantitative analysis of intra- and inter-specific 


been used to evaluate cranial, 
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variations at the morphological level is useful for 


detecting patterns of geographic variations and 
delimiting intra- or inter-specific evolutionary units. To 
date, however, there are currently no published reports of 
quantitative analysis based on morphological characteristics 
that would allow the identification of the morphotypes in 
the complex of P. philippensis and P. petaurista. 

To discuss the taxonomic relationships of P. 
and P. 


petaurista and test previous taxonomic hypotheses, the 


philippensis, P. yunanensis, P.  hainana, 
present study conducted a comprehensive morphometric 
study on the above Chinese Petaurista species based on 
samples subsequently collected from southwest Yunnan 
and the Island of Hainan, China. Multivariate analyses 
were used to produce an overview of the associations 
between morphological variables and species patterns 
and discuss the taxonomic implications of these flying 
morphometric study could be 

studies of variations of DNA 
sequences in flying squirrels. 


1 Materials and Methods 


squirrels. Our 
complementary to 


1.1 Specimens and data collection 

According to the taxonomic assignments of Allen 
(1940) and Zhang et al (1997), a total of 60 intact adult 
skull specimens of P. petaurista, P. yunanensis, P. 
hainana, and P. philippensis were examined for 
morphometric study (Append. I). These specimens are 
from the Kunming Institute of Zoology, Chinese 
Academy of Sciences (KIZ, CAS) (Kunming, China), the 
Institute of Zoology, Chinese Academy of Sciences (IOZ, 
CAS) (Beijing, China), and the Guangdong 
Entomological Institute (GDEI) (Guangzhou, China). 

Twenty-six cranial variables taken with a digital 
caliper to the nearest 0.01 mm were used in the 
morphometric analysis as described by Musser (1979), 
Musser & Heaney (1992), Xia et al (2006), and Yang et 
al (2005), and depicted in Fig. 1 following Huang’s 
description (1995). The variables measured included: 
maximum length of skull (GLS), condylobasal length 
(CBL), basal length (BL), occipito-nasal length (ONL), 
palatal length (PL), length of palatal bridge (PBL), 
length of upper tooth row (LUTR), length of upper 
molars (LUM), maximum upper molars breadth 
(GUMB), rostral length (ROL) and breadth (ROB), 
auditory bulla length (ABL) and breadth (ABB), breadth 
of zygomatic plate (BZP), breadth of occipital condyles 
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Fig. 1 Illustration of the 26 cranial variables used in the study 


(BOO), height of occipital (HO), zygomatic breadth 
(ZOB), mastoid breadth (MTB), nasal length (NL) and 
breadth (BN), mandible length (ML), height of mandible 
(THM), length of lower diastema (LLD), length of lower 
molar row (LLMR), length of lower tooth row (LTR), 
and mandibular height (MH). In addition, the head and 
body length (HB), tail length (TL), hind foot length 
(HFL), and ear length (EL), which were compared to the 
original measurements labeled on the skins by the 
collectors. 
1.2 Data analysis 

Statistical analyses were performed using SPSS 
version 11.0 (SPSS Inc., Chicago, IL, USA). All 
variables were transformed into logarithms to eliminate 
the bias effect of large measurements in multivariate 
(D’Elia & Pardiñas, 2004). Statistical 
differences were considered significant at P<0.05. 


analysis 


In this study, all related data were subjected to one- 





way ANOVA for calculating mean+SD. T-test was used 
to assess the sexual dimorphism between male and 
female groups by comparing the group means of cranial 
measurements. Multiple comparisons between taxa were 
made for all 26 cranial measurements to evaluate 
variations between samples. Multivariate analyses, 


including principal components analysis (PCA) and 
discriminant function analysis (DFA), were carried out to 
evaluate the degree of similarity and dissimilarity in 
cranial structures between the putative species and to 
determine how the taxa were related when all cranial 
characters measured are considered simultaneously. 

The PCA is based upon the variance-covariance 
matrix of the log-transformed variables. The eigenvector 
scores describing the relative significance of each 
variable to the principal components were used to 
compare the cranial morphological similarities and 
differences. The PCA scatter-plot visually represented 
the variation among different individuals of the samples. 
The DFA was performed to investigate the integrity of 
the pre-defined groups and to predict group membership 
of specimens with the linear models of variables. Based 
on the derived discriminant functions, each individual 
was allocated to the group with nearest centroid, and the 
proportion of individuals allocated to each group was 
calculated. 


2 Results 





Mean+SD of 4 external and 26 cranial variables for 
the four taxa are presented in Tab. 1. 
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Tab. 1 Four external and 26 cranial measurements of four species of Petaurista: (mean+SD)/range 
GLS CBL ONL PL PBL LUTR LUM 
P. yunanensis 771.16+2.73 71.55+2.51 66.29+2.58 76.9042.70 40.08+1.72 26.97+1.06 37.57+£1.36 12.4640.42 
n=15 71.67-81.56 67.79-75.68 62.24-70.53 71.64-81.09 37.55-43.05 25.05-28.59 35.73-39.80 11.68-13.06 
P. philippensis 79.10+2.01 713.93+2.14 68.60+2.09 78.70+2.08 41.82+1.47 28.10+1.32 38.94+1.13 12.46+0.47 
n=21 75.35-82.71 69.10-77.75 64.44-72.39 74.96-82.77 39.11-43.80 26.13-30.64 37.27-40.77 11.67-13.29 
P. hainana 713.45+1.90 68.92+1.80 63.83+1.77 72.82+1.83 39.05+1.20 25.35+0.92 36.14+1.12 11.16+0.31 
n=17 68.96-76.68 65.19-72.04 60.34-66.86 68.59-75.52 36.93-41.25 23.96-27.29 34.04-37.99 10.58-11.78 
P. petaurista 63.2542.12 59.40+2.00 55.05+1.84 62.8442.24 33.64+1.54 22.23+1.25 30.88+0.98 9.92+0.43 
n=7 60.16-66.43 56.98-62.67 52.69-58.30 59.39-66.14 32.26-35.86 20.25-23.80 29.72-32.34 9.15-10.42 
GUMB ROL ABL ABB BZP BOO HO 
P. yunanensis 20.25+0.89 16.61+1.07 15.64+0.82 13.2840.72 10.17+0.45 8.20+0.66 18.2740.77 8.43+0.79 
n=15 18.70-21.46 14.94-18.14 13.94-16.85 12.22-14.62 8.83—10.61 6.94-9.12 17.00-19.67 6.87—9.76 
P. philippensis 20.43+0.79 17.63+0.81 16.56+0.72 14.81+0.67 11.18+0.48 7.73+0.43 18.58+0.88 8.53+0.57 
n=21 19.07-22.57 16.05-18.93 15.36-18.11 13.73-15.76 10.50-12.41 6.60-8.74 17.48-20.22 7.47-9.57 
P. hainana 19.26+0.72 15.88+0.76 16.12+0.73 14.27+0.36 10.53+0.63 6.67+0.46 17.48+0.82 9.45+0.69 
n=17 18.14-20.57 14.06-17.34 14.51-17.30 13.60-14.91 9.36-11.48 5.96-7.66 15.57-18.62 8.27-10.72 
P. petaurista 16.60+40.46 12.72+0.51 13.81+0.65 13.1940.38 9.33+0.48 6.04+0.35 14.39+0.53 7.61+0.24 
n=7 16.05-17.46 12.09-13.30 12.65-14.53 12.81-13.88 8.54-9.78 5.53-6.53 13.51-15.00 7.34-8.04 
ZOB MTB BN ML THM LLD LLMR 
P. yunanensis 51.00+1.74 29.34+1.11 23.82+1.56 15.15+1.00 53.96+2.09 31.88+1.19 10.30+0.81 13.66+0.45 
n=15 48.18-54.40 27.26-31.24 20.59-26.05 13.94-16.91 49.54-57.31 29.95-34.00 8.95-12.01 13.13-14.66 
P. philippensis 51.424+1.72 30.05+1.31 24.34+1.46 15.30+40.96 54.82+1.33 32.23+1.23 10.82+0.65 14.03+0.56 
n=21 46.48-54.73 27.71-32.11 21.98-27.34 13.53-17.40 52.31-57.68 29.61-34.64 9.42-11.94 13:10-15:37 
P. hainana 49.09+1.24 27.58+1.07 22.76+1.23 14.61+0.61 50.92+1.21 30.08+1.25 10.21+0.53 12.53+0.35 
n=17 46.79-51.34 25.96-29.39 20.89-24.71 13.38-15.50 48.81-52.71 27.96-32.40 8.97—10.88 11.83-13.12 
P. petaurista 41.34+1.10 22.67+0.92 19.47+1.31 11.13+40.47 43.42+1.11 26.08+1.07 8.52+0.42 10.74+0.26 
n=7 39.16—42.50 21.54-23.65 17.84-21.15 10.39-11.71 42.08-44.93 24.65-28.04 7.82-9.20 10.39-11.16 
LTR MH HB TL HFL EL 
K 443.07+26.80 528.69+49.19 82.93+4.91 44.14+3.32 
Feyunanensiz EL 12 21EG 420.00-520.00 440.00-620.00 75.00-92.00 40.00-52.00 
n=15 33.95-37.19 11.38-13.80 a i 5 i 
n=14 n=13 n=14 n=14 
D f 459.53+33.17 549.95+41.89 83.42£7.74 49.21+4.20 
P. philippensis 36-93£1:18 12112032 384.00-530.00 470.00—620.00 65.00—100.00 42.00-58.00 
n=21 34.98-39.89 11.02-12.95 
n=19 n=19 n=19 n=19 
é 413.63427.29 486.31430.60 77.004£5.21 41.59+4.12 
P. hainana 34.2140.87 11.61+0.41 
n=17 32.35-35.41 10.88-12.32 348.00=450:00 S000; 320,00 Sicha OOO ADI 
n=16 n=16 n=17 n=17 
337.50+3.54 360.00+0.00 60.00+0.00 42.00+0.00 
P. petaurista 28.96+0.83 9.77+0.25 A a À D 
n=7 27.81-29.90 9.35-10.06 BAA tad POR T ERIN ua papas 


Variable codes are given in the text and Fig. 1. 


2.1 Univariate analysis 

Univariate comparison revealed that the means of 
all variables were significantly different and, in general, 
tended to become progressively larger from P. petaurista, 
P. hainana, P. yunanensis, to P. philippensis. The t-tests 
of Equality of Group Means on 54 (30 males, 24 females) 
out of 60 specimens indicated there was no sexual 
dimorphism in the 26 cranial variables in the four 
Petaurista groups (Tab. 2). Quantitative pairwise 
comparisons of all cranial variables between taxa 
indicated that P. yunanensis was morphologically similar 
to P. philippensis, with 11 cranial measurements 
showing no significant difference (P>0.05) (Tab. 3). 
Also, five cranial variables were not significantly 
different between P. yunanensis and P. hainana. 


2.2 Multivariate analysis 

In PCA, the eigenvalues for the first three principal 
components were 20.62, 1.64 and 0.85, respectively, 
accounting for 88.91% of the total variance (Tab. 4). 
Most characteristics with high positive loadings on the 
first principal component suggested that this component 
(79.32% of the total variance) represented size variation 
within the samples. All specimens on the first principal 
component were clustered as three groups, P. petaurista, 
P. hainana, and the group of P. yunanensis and P. 
philippensis with considerable overlaps. The second 
principal component (6.32% of the total variance) was 
strongly correlated with ROB, ABL, ABB, and BZP 
(loadings>0.50), and the third principal component 
(3.27% of the total variance) was correlated primarily 
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Tab. 2 t-tests of Equality of Group Means for male and 
female (variable codes are given in the text and Fig.1) 





with HO (loadings>0.50) (Tab. 4). The first two principal 
components separated all specimens as four distinct 








Wilks' Lambda F dfl df2 P groups (Fig. 2). 
S ees ees ; 2 ae The DFA identified the major patterns of 
CBL 0.976 1.274 1 52 0.264 
BL 0.980 1.087 1 52 0.302 morphological divergence in the crania among the four 
ONL 0.981 1.021 1 52 0.317 Petaurista groups. The variation pattern reflected by the 
ae mee en F i first two discriminant functions was consistent with the 
LUTR 0.993 0.386 1 52 0.537 morphological variations observed in the PCA, and all 
LUM 9224 0:237 i 32 0-588 samples were clearly clustered as four distinguishable 
GUMB 0.983 0.924 1 52 0.341 a Hone) S : 
ROL 1.000 0.006 1 52 0.939 groups based on the 1° and 2™ discriminant functions 
ROB 0.959 2.214 1 52 0.143 (Fig. 3). In specimen reclassification by DFA, all 
ABE 0286 0:746 l 2 9.322 individuals were properly assigned to their original 
ABB 0.978 1.191 1 52 0.280 
BZP 0.956 2.368 1 52 0.130 groups on the basis of the studied measurements. Fig. 4 
BOO 0.993 0.377 l 52 0.542 is the geographic distributions of all Petaurista samples 
HO 0.996 0.233 1 52 0.632 SAN: 
ZOB 0.967 1.775 1 52 0.189 used in this study. 
MTB 0.993 0.359 1 52 0.552 x , 
NL 0.998 0.106 1 52 0.746 3 Discussion 
BN 0.994 0.291 1 52 0.592 
ML 0.992 432 1 52 0.514 One contentious issue regarding the Chinese 
THM 0.955 2.438 1 52 0.124 Petaurista is the taxonomic status of P. yunanensis, P. 
LLD 0.986 0.739 1 52 0.394 philippensis, P. hainana, and the populations of P. 
ry aes UKE ' a pass petaurista in China, which have long been controversial 
MH 0.981 0.988 1 52 0.325 (Corbet & Hill, 1992; Ellerman, 1940; Ellerman & 
P<0.05 Morrison-Scott, 1950; Hoffmann et al, 1993; Huang et al, 
Tab. 3 Multiple Comparisons between the four study species 
P. yunanensis vs P. yunanensis vs P. yunanensis vs P. philippensis vs P. philippensis vs P. hainana vs 
Variables P. philippensis P. hainana P. petaurista P. hainana P. petaurista P. petaurista 
Mean Mean Mean Mean Mean Mean 
Difference Difference Difference Difference Difference Difference 
GLS 0.011* 0.015 0.021** 0.000 ~—:0.086** 0.000  0.032** 0.000 0.097** 0.000  0.065** 0.000 
CBL 0.016* 0.002 0.016* 0.003  0.081** 0.000  0.033** 0.000  0.097** 0.000  0.065** 0.000 
BL 0.015* 0.003 0.016* 0.002  0.081** 0.000  0.031** 0.000  0.095** 0.000  0.064** 0.000 
ONL 0.010* 0.027 0.024** 0.000  0.088** 0.000  0.034** 0.000  0.098** 0.000  0.064** 0.000 
PL 0.019* 0.001 0.011 0.058  0.076** 0.000  0.030** 0.000 0.095** 0.000  0.065** 0.000 
PBL 0.018* 0.008  0.027** 0.000  0.084** 0.000  0.045** 0.000  0.102** 0.000  0.057** 0.000 
LUTR 0.015* 0.002 0.017* 0.001  0.085** 0.000  0.032** 0.000 0.100** 0.000 0.068** 0.000 
LUM 0.001 0.911 0.048** 0.000  0.099** 0.000  0.047** 0.000 0.099** 0.000 0.051** 0.000 
GUMB 0.003 0.542 0.022* 0.001  0.086** 0.000  0.025** 0.000 0.089** 0.000  0.064** 0.000 
ROL 0.026* 0.001 0.019* 0.018 0.116%" 0.000  0.046** 0.000  0.142** 0.000  0.096** 0.000 
ROB 0.025* 0.001 0.013 0.070  0.054** 0.000 0.011 0.103 0.078** 0.000  0.067** 0.000 
ABL 0.047** 0.000  0.032** 0.000 0.002 0.771 0.015* 0.012 0.050** 0.000  0.034** 0.000 
ABB 0.041** 0.000 0.014 0.065  0.038** 0.000  0.027** 0.000  0.079** 0.000  0.052** 0.000 
BZP 0.025* 0.015 0.089** 0.000  0.132** 0.000  0.065** 0.000  0.107** 0.000 0.043* 0.002 
BOO 0.007 0.293 0.019* 0.007  0.104** 0.000  0.027** 0.000  0.111** 0.000  0.084** 0.000 
HO 0.006 0.578 0.051** 0.000 0.043* 0.005  0.044** 0.000 0.049* 0.001  0.094** 0.000 
ZOB 0.003 0.462 0.016* 0.001  0.091** 0.000  0.020** 0.000  0.094** 0.000  0.075** 0.000 
MTB 0.010 0.090  0.027** 0.000  0.112** 0.000  0.037** 0.000 0.122** 0.000  0.085** 0.000 
NL 0.009 0.311 0.019* 0.044  0.088** 0.000 0.029* 0.002  0.097** 0.000  0.068** 0.000 
BN 0.004 0.600 0.015 0.083 0.133** 0.000 0.020* 0.017  0.138** 0.000 0.118** 0.000 
ML 0.007 0.114  0.025** 0.000  0.094** 0.000  0.032** 0.000 0.101** 0.000  0.069** 0.000 
THM 0.005 0.383 0.025** 0.000  0.087** 0.000  0.030** 0.000 0.092** 0.000 0.062** 0.000 
LLD 0.022* 0.023 0.003 0.767 0.082** 0.000 0.025* 0.008  0.103** 0.000  0.079** 0.000 
LLMR 0.012* 0.018  0.037** 0.000  0.104** 0.000  0.049** 0.000  0.116** 0.000  0.067** 0.000 
LTR 0.016* 0.001 0.018** 0.000  0.090** 0.000  0.034** 0.000 0.106** 0.000  0.072** 0.000 
MH 0.003 0.583 0.022* 0.002 0.096** 0.000  0.025** 0.000  0.100** 0.000 0.075** 0.000 


*: P <0.05; **: P <0.001 


124 Zoological Research Vol. 33 


Tab. 4 Factor loadings and percentage of variance 


















































explained for principal component analysis Q pi 
(variable codes are given in the text and Fig. 1) A O & 
Variables PC1 PC2 PC3 Pp O 
GLS 0.898 0.347 0.221 K & 
CBL 0.868 0.379 0.212 a A 
N 
BL 0.876 0.388 0.227 9 Aa ê a 
A 
ONL 0.904 0.336 0.201 a 
PL, 0.841 0.422 0.224 D P. yunanensis o 
(m 
PBL 0.872 0.347 0.024 © P. philippensis 
LUTR 0.889 0.375 0.192 4 P. hainana 
LUM 0.944 0.164 -0.016 + P. petaurista 
MB 0.839 0.250 0.330 -3 : : 5 : 
or -3 -2 -1 0 1 2 
ROL 0.820 0.420 0.151 PCI 
ROB 0.535 0.586 0.470 ; i 
ABL 0.129 0.895 0.207 Fig. 2 Scatterplots of the samples in PCA space 
ABB 0.393 0.766 0.085 6 
BZP 0.833 -0.608 -0.028 
BOO 0.847 0.279 0.258 4 mm 
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Fig. 4 Geographic distributions of samples used in the study 
A: Dulongjiang; B: Bijiang; C: Yunlong; D: Yingjiang; E: Lianghe; F: Tengchong; G: Longling; H: Luxi; I: Cangyuan; J: Xishuangbanna; K: Xinping; L: Mile; 
M: Lvchun; N: Pingbian; O: Xiangzhou; P: Northern Guangdong; Q: Hainan. 
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1995; Oshida et al, 2000a, 2000b; Thorington & 
Hoffmann, 2005; Wang, 2003; Yu et al, 2006). By using 
a considerably extended set of morphometrical 
applying 
multivariate morphometric analyses, results of the 


characters (26 cranial variables) and 
present study confirmed the significant craniological 
differences in P. petaurista, P. hainana, P. yunanensis, 
and P. philippensis, with P. petaurista having the most 
pronounced morphological variations, particularly in 
metrical components of cranial and body size. 

Pelage coloration had been applied for classification 
of flying squirrels and led to many taxonomical 
disagreements due to numerous color variations within 
Petaurista, even between different sexes (Allen, 1940; 
Ellerman & Morrison-Scott, 1950; Oshida et al, 2004b). 
A series of color variations in pelage were observed 
among the Petaurista forms including both sexes (Allen, 
1940; Oshida et al, 2004a), but the pairwise comparison 
of each of the 26 cranial variables revealed no sexual 
dimorphism in any of the four groups, implying that the 
divergence of coloration patterns in forms was due to 
environmental or genetic fluctuations over time. 

Both univariate and multivariate analysis revealed 
that the skull morphometric characters used in this study 
were effective for discriminating the four Petaurista 
groups. Our analyses demonstrated that the four 
morphotypes of Petaurista were distinguished by a 
number of cranial characteristics. The morphometric 
variables, which caused the major distinction between 
those groups, were specifically located in the occipital, 
supraocular, and rostral regions, as well as in the prootic- 
squamosal length. Patterns of molecular sequence 
variations from previous studies (Oshida et al, 2000a; Yu 
et al, 2006) and the cranial morphological differences 
observed in this study indicate that P. hainana, P. 
philippensis, and P. yunanensis could be recognized as 
three distinct species. These differences were clear and 
reinforced the existence of three morphotypes of the 
complex P. philippensis; although the degree and form of 
the morphological differences might be related to their 
geographical variations. Even though molecular data 
suggested that P. philippensis is closely related to P. 
hainana and significantly distinct from P. yunanensis 
(Yu et al, 2006), many characteristics beyond those 
related to external morphology and pelage coloration 
were observed to be held in common between P. 
yunanensis and P. philippensis. A good example is the 
similarity in the pattern of overall cranial structure in 


multiple comparison analysis. The degree of the 
morphological variations between P. yunanensis and P. 
philippensis was less than that observed between P. 
philippensis and P. hainana, with eleven cranial 
measurements showing no significant difference (P>0.05) 
(Tab. 3). Petaurista yunanensis occurs from extreme 
southwestern Yunnan into Myanmar and Indochina and 
is extensively sympatric with P. philippensis in 
southwestern China (Wang, 2003; Zhang et al, 1997). 
The sharing of morphological characteristics between P. 
philippensis and P. yunanensis is related to their similar 
living conditions. 

Petaurista hainana was considered a valid species 
based on both molecular and morphological data (Huang 
et al, 1995; Yu et al, 2006). Our morphometric results 
were concordant with mtDNA data of previous research 
(Oshida et al, 2000a; Yu et al, 2006) and demonstrated 
the significant differences between P. hainana and P. 
yunnanensis/P. philippensis. In both PCA and DFA, P. 
hainana was clearly separated from other three groups 
(Fig. 2, 3), with 21/26 cranial variables being 
significantly different (P<0.05) (Tab. 3). Petaurista 
hainana is confined to tropical forests on Hainan Island 
of China and P. philippensis and P. yunanensis are 
distributed in mountainous coniferous, dry deciduous and 
evergreen forests at different elevations in western 
Yunnan of China. The phenotypic divergence of P. 
hainana in relation to P. philippensis and P. yunanensis 
is likely associated with their geographical distributions 
and living conditions and could be viewed as a reflection 
of adaptations to various ecological niches. The 
differences in skull morphology suggest that P. hainana 
is neither the synonym of P. philippensis nor a 
subspecies of P. yunanensis or P. petaurista (Corbet & 
Hill, 1992; Thorington & Hoffmann, 2005; Wang, 2003), 
but a valid species in its own right. 

The greatest distinction observed was between P. 
petaurista and other Petaurista forms. Petaurista 
petaurista displayed a relatively high level of diversity in 
with 22/26 
significantly different from P. hainana, P. philippensis, 
and P. yunanensis at P<0.001 level (Tab. 3). Based on 26 
morphological cranial variables, the specimens of P. 


skull morphology, cranial variables 


petaurista formed a distinct aggregate in both PCA and 
DFA (Fig. 2,3), consistent with the results of Oshida et al. 
(2000a) and Yu et al (2006). It is obvious that P. 
petaurista, P. hainana, P. philippensis, and P. yunanensis 
are taxonomically distinct and distinct valid species. 
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Append I: specimens examined 

Petaurista yunanensis n=15 

Yingjiang, Yunnan: IOZ 25849(3). Gongshan, 
Yunnan: KIZ 73442(2), 73445(3), 73744(3), 73745(9), 
73823(3), 830207(), 90039(9), 90043(2), 90051(8), 
90407. Tengchong, Yunnan: KIZ 76348(@). Lianghe, 
Yunnan: KIZ 650236(4), 650237(3). Bijiang, Yunnan: 
KIZ 780102(). 

Petaurista philippensis n=21 

Xishuangbanna, Yunnan: IOZ 10457(8), 10458(8), 
10460(2), 15041(Q), 15042(2), 


References: 


Allen GM. 1940. The mammals of China and Mongolia[M]. New York: 
The American Museum of Natural History, 729-745. 

Corbet GB, Hill JE. 1992. The Mammals of the Indomalayan Region: A 
Systematic Review[M]. Oxford: Oxford University Press, 308-313. 

D’Elia G, Pardiñas UFJ. 2004. Systematics of Argentinean, Paraguayan, 
and Uruguayan swamp rats of the genus Scapteromys (Rodentia, 
Cricetidae, Sigmodontinae)[J]. J. Mammal, 85(5): 897-910. 

Ellerman JR. 1940. The Families and Genera of Living Rodents. 
Volume I[M]. London: British Museum, 281-290. 

Ellerman JR, Morrison-Scott TCS. 1950. Checklist of Palaearctic and 
Indian mammals, 1758 to 1946[M]. London: British Museum 
(Natural History), 460-465. 

Hoffmann RS, Anderson CG, Thorington RW, Jr., Heaney LR. 1993. 
Family sciuridae[M]. //Wilson DE, Reeder DM. Mammal Species 
of the World: A Taxonomic and Geographic Reference. 2" ed. 
Washington D.C.: Smithsonian Institution Press, 462-463. 

Huang WJ, Chen YX, Wen YX. 1995. Rodentia of China[M]. Shanghai: 
Fudan University Press, 50-59. (in Chinese) 

Muñoz-Muñoz F, Perpiñán D. 2010. Measurement error in morphometric 
studies: comparison between manual and computerized methods[J]. 
Ann Zool Fennici, 47: 46-56. 

Musser GG. 1979. Results of the Archbold Expeditions. No.102. The 
species of Chiropodomys, arboreal mice of Indochina and the 
Malay Archipelago[J]. Bull Am Nus Nat Hist, 162: 377-445. 

Musser GG, Heaney LR. 1992. Philippine Rodents: Definitions of 
Tarsomys and Limnomys plus a preliminary assessment of 
phylogenetic patterns among native Philippine murines (Murinae, 
Muridae)[J]. Bull Am Nus Nat Hist, 211: 1-138. 

Oshida T, Lin LK, Masuda R, Yoshida MC. 2000a. Phylogenetic 
relationships among Asian species of Petaurista (Rodentia, 
Sciuridae), inferred from mitochondrial cytochrome b gene 
sequences [J]. Zool Sci, 17(1): 123-128. 

Oshida T, Lin LK, Yanagawa H, Endo H, Masuda R. 2000b. 
Phylogenetic relationships among six flying squirrel genera, 
inferred from mitochondrial cytochrome b gene sequences[J]. Zool 
Sci, 17(4): 485-489. 

Oshida T, Shafique CM, Barkati S, Fujita Y, Lin LK, Masuda R. 2004a. 


Vol. 33 


15043(3), 15044(3), 61-003(2). Luxi, Yunnan: 
IOZ 24009(2), 24010(9). Mile, 

Yunnan: KIZ 84255(3), 84258(3), 84300(<). 
Lvchun, Yunnan: KIZ 72104(). 

Longling, Yunnan: KIZ 620028(¢). Cangyuan, 
Yunnan: KIZ 78053(<). Pingbian, 

Yunnan: KIZ 84117(3), 
Yunnan: KIZ 77017(Q). Yunlong, 

Yunnan: KIZ 200369(9). Yunnan: KIZ 92005. 

Petaurista hainana n=\7 

Hainan: GDEI 0403(2), 0404(<), 
0524(3), 0611(), 0618(4), 0621(), 

0622(3), 0623(2), 0624(4), 0625(3), 0626(<4), 
0703(3), 0704(3), 0705(), 

0714(8), 0717(2). 

Petaurista petaurista n=7 

Xiangzhou, Guangxi: IOZ 00433(4). Northern 
Guangdong: GDEI 2279, 2284, 2286(Q), 2287, 2288(3), 
2289. 


84119(4). Xinping, 


0499(9), 


A preliminary study on molecular phylogeny of giant flying 
squirrels, genus Petaurista (Rodentia, Sciuridae) based on 
mitochondrial cytochrome b gene sequences[J]. Russian Journal 
of Theriology, 3: 15-24. 

Oshida T, Shafique CM, Barkati S, Yasuda M, Hussein NA, Endo H, 
Yanagawa H, Masuda R. 2004b. Phylogenetic position of the 
small Kashmir flying squirrel, Hylopetes fimbriatus (= 
Eoglaucomys fimbriatus), in the subfamily Pteromyinae[J]. Can J 
Zool, 82(8): 1336-1342. 

Slabova M, Frynta D. 2007. Morphometric variation in nearly 
unstudied populations of the most studied mammal: the non- 
commensal house mouse (Mus musculus domesticus) in the near 
East and Northern Africa[J]. Zool Anz, 246(2): 91-101. 

Thorington RW Jr, Hoffmann RS. 2005. Family Sciuridae[M] // Wilson 
DE, Reeder DM. Mammal Species of the World: A Taxonomic and 
Geographic Reference. 3™ ed. Washington D.C.: The Johns 
Hopkins University Press, 770-772. 

Wang YX. 2003. A Complete Checklist of Mammal Species and 
Subspecies in China: A Taxonomic and Geographic Reference[M]. 
Beijing: China Forestry Publishing House, 155-159. (in Chinese) 

Xia L, Yang QS, Ma Y, Feng ZJ, Zhou LZ. 2006. Guide to the 
measurement of mammal skull III: Rodentia and Lagomorpha{J]. 
Chinese J Zool, 41(5): 68 -71. (in Chinese) 

Yang QS, Xia L, Ma Y, Feng ZJ, Quan GQ. 2005. A guide to the 
measurement of mammal skull I: Basic measurement[J]. Chinese J 
Zool, 40(3): 50-56. (in Chinese) 

Yu FR, Yu FH, Peng JF, Kilpatrick CW, McGuire PM, Wang YX, Lu 
SQ, Woods CA. 2006. Phylogeny and biogeography of the 
Petaurista philippensis complex (Rodentia: Sciuridae), inter- and 
intraspecific relationships inferred from molecular and 
morphometric analysis [J]. Mol Phylogenet Evol, 38(3): 755-766. 

Zelditch ML, Swiderski DL, Sheets HD, Fink WL. 2004. Geometric 
Morphometrics for Biologists: A Primer[M]. New York: Elsevier 
Academic Press. 

Zhang YZ et al. 1997. Distribution of Mammalian Species in China[M]. 
Beijing, China: China Forestry Publishing House, 150-154. (in 
Chinese) 


